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First Calculation of Hyperon Axial Couplings from Lattice QCD 
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In this work, we report the first lattice calculation of hyperon axial couplings, using the 2-|-l-flavor 
MILC configurations and domain-wall fermion valence quarks. Both the E and H axial couplings 
are computed for the first time in lattice QCD. In particular we find that = 0.450(21)stat (27)syst 
and SEE = -0.277(I5)stat (19),yst. 
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I. INTRODUCTION 

In recent years lattice QCD calculations have demon- 
strated remarkable progress in computing hadron proper- 
ties from first principles. In particular, hadron structure 
has been a major focus of the lattice QCD community. 
Recently the nuclcon axial coupling gA has been com- 
puted [J, i with good precision and has been shown to 
be in agreement with the very well known experimental 
result within the systematic and statistical errors of the 
calculation. The success of such calculations motivates 
us to study the axial couplings for all the octet baryons. 
The experimental knowledge of such couplings is not as 
good as in the case of gA and in certain cases, such as the 
axial couplings gss and ^hh, their values are not known 
experimentally, and theoretical estimates are rather im- 
precise. 

The hyperon axial couplings arc important parameters 
entering the low-energy effective field theory description 
of the octet baryons. At the leading order of SU(3) heavy 
baryon chiral perturbation theory, these coupling con- 
stants are linear combinations of the universal coupling 
constants D and F, which enter the chiral expansion of 
every baryonic quantity, including masses and scattering 
lengths. These coupling constants are needed in the effec- 
tive field theory description of both the non-leptonic de- 
cays of hyperons, and the hyperon- nucleon and hyperon- 
hyperon scattering phase shifts Hyperon- nucleon 

and hyperon-hyperon interactions are essential in under- 
standing the physics of neutron stars where hyperon and 
kaon production may soften the equation of state of dense 
hadronic matter. 

Studying the octet baryon axial couplings on the lat- 
tice is no more complicated than computing the nuclcon 
axial coupling gA- The lack of experimental informa- 
tion in cases such as gss and gss gives us the opportu- 
nity to make predictions using lattice QCD. Previously, 
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there have been two attempts to determine these coupling 
constants using different theoretical approaches: chiral 
perturbation theory H and the large- TVc limit 0]. 
M. J. Savage ct al. 0] use chiral perturbation theory 
to work out the one-loop corrections due to SU(3) sym- 
metry breaking, and predict 0.35 < gss < 0.55. Simi- 
larly, 0.18 < -gEE < 0.36, after taking the SU(3) limit, 
5es = F, which is 0.5 at the leading order. Using the 
large- iVj;. approach, Dai et al. Q and Flores-Mendieta 
et al. [y| predict a range of 0.30 < gss < 0.36 and 
0.26 < —gEE < 0.30. Both approaches give very loose 
bounds on the values of these coupling constants; hence, 
a lattice QCD calculation has the opportunity to make a 
substantial improvement. 



Lattice QCD calculations can now provide much more 
stringent theoretical estimates of these axial couplings. 
The various systematic errors that enter such calcula- 
tions can be controlled, giving us the ability to compute 
the predictions of QCD very precisely. The systematic 
errors that we need to control are due to the finite lat- 
tice spacing, the finite volume and chiral extrapolations. 
Lattice calculations arc performed on a discrete space- 
time in a finite volume, using Monte Carlo integration to 
directly evaluate the path integral; continuum physics is 
recovered by taking the lattice spacing to zero (a — > 0) 
and the volume to infinity (V oo). In addition, using 
current computer resources, we cannot yet calculate at 
the physical pion mass. Using chiral perturbation theory 
and calculations at multiple heavier pion masses which 
are more affordable, we can extrapolate quantities of in- 
terest to the physical limit. Such calculations also help to 
determine the low-energy constants of the chiral effective 
theory and allow us to study the quark-mass dependence 
of our observables. 

The structure of this article is as follows: In Scc.|TTl we 
describe our lattice setup, operators and the parameters 
of the calculation. In Sec. IIIIl we discuss the chiral ex- 
trapolation, and finally, our conclusions are presented in 
Sec. lIVl Our preliminary results can be found in Ref. 0]. 
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II. THE LATTICE CALCULATION 

The axial coupling constants arc defined as the zero 
momentum transfer limit of the axial form factor Ga(9^) 
that parametratizes the matrix element (B\0\B) where 
B is one of the baryons N, S, S, and O is the local axial 
current (97^759)- This matrix element takes the form 



{B\A^{q)\B) = ub{p') 



2Mf 



xuB{p)e-''>-\ 



(1) 





mOlO 


m020 


m030 


m040 


m050 


(MeV) 


354.2(8) 


493.6(6) 


594.2(8) 


685.4(19) 


754.3(16) 




2.316(7) 


3.035(7) 


3.478(8) 


3.822(23) 


4.136(20) 


TTIk/ fiT 


3.951(14) 


3.969(10) 


4.018(11) 


4.060(26) 


4.107(21) 


confs 


612 


345 


561 


320 


342 


9A,N 


1.22(8) 


1.21(5) 


1.195(17) 


1.150(17) 


1.167(11) 


ffSE 


0.418(23) 


0.450(15) 


0.451(7) 


0.444(8) 


0.453(5) 




-0.262(13) 


-0.270(10) 


-0.269(7) 


-0.257(9) 


-0.261(7) 



TABLE I: Results of our calculation 



where is the Dirac spinor. Gp is the induced pseu- 
doscalar form factor and q is the momentum transfer. 
The axial charge for baryon B is defined as qa bb = 
GA[q^ = 0). 

On the lattice, we can extract the matrix element 
{B by calculating a three-point function using 

zero initial and final momentum for the baryon states. 
This matrix element needs to be rcnormalized, because 
we use the local axial current in order to simplify our cal- 
culation. The renormalization constant Za can be eas- 
ily computed through two-point meson correlation func- 
tions. Details of our lattice formulation and methods can 
be found in Refs. 

For the rest of the paper, we will define the rcnor- 
malized zero mometum transfer matrix elements as fol- 
lowing: for the nucleon qa = Za{N\A^\NY''^, the E 
3es = (E l^^ l S])''^V2 (the factor of 2 is coming from a 
Clebsch-Gordan coefficient so that gss = F in the SU(3) 
limit) and for the S ^hh = Za{^ S)''^*. 

In this calculation, we use (improved) staggered 
fermion action (asqtad) [l^, [HI [l^ for the sea quarks 
and domain- wall fermions (DWF) [13, [ll [H, [il| for the 
valence sector. This way we take advantage of the pub- 
licly available gauge configurations generated by MILC 
collaboration, only having to compute the valence quark 
propagators needed for the matrix elements. In ad- 
dition, domain-wall fermions used in the valence sec- 
tor are automatically 0{a) improved and have chiral 
and flavor symmetry which simplifies operator mixing, 
renormalization, and chiral extrapolation at finite lattice 
spacing [13, [H, [H, [l^l, making them particularly well 
suited for the purpose of our calculation. This mixed ac- 
tion approach has been successfully employed by LHPC 
and NPLQCD for c omp utations of nucleon matrix el- 
ements [U, [23, [23l. l24i as well as scatte ring lengths, 
decay constants fil,!!!, [13, [M [H, [M, [M [12 and 
masses [33^^, ,35|] ^. 



Staggered fermions come in four "tastes" which must be removed 
by taking fractional powers of the fermionic determinant. Al- 
though no theoretical proof of the validity of this approach ex- 
ists, there is significant evidence in the literature that supports 
the conjecture that the continuum limit of such formulation is 
QCD. For a recent review see Ref. [ssll and references therein. 



The gauge configurations are generated with 2-1-1 fla- 
vors of staggered fermions (configuration ensembles gen- 
crated by the MILC collaboration [13] )■ The pion mass 
ranges from 350 to 750 MeV in a lattice box of size 2.6 fm. 
The gauge fields that enter the domain-wall fermion ac- 
tion are hypercubic smeared to improve the chiral sym- 
metry, and gauge invariant Gaussian smearing has been 
used for the interpolating operators to improve the signal 
(see [21I, [H, [13] for details) . The source-sink separation 
is fixed at 10 time units. The number of configurations 
used from each ensemble ranges from 300 to 600. Our 
results are presented in Table IH 



III. HYPERON AXIAL COUPLING 
CONSTANTS 

A. SU(3) symmetry breaking 

One way to probe SU(3) symmetry breaking in the 
axial couplings is to monitor the quantity 5su(3)! defined 
as 

<5su(3) = 5A - 2.0 X gss + .gsE = ^ c„x"; (2) 



where x is (m^ — m^)/(47r/^). Figure [T] shows ^su(3) 
a function of x. Note that the value increases monoton- 
ically as we go to lighter pion masses. Our lattice data 
suggest that a (5su(3) ~ dependence is strongly pre- 
ferred, as the plot of (5su(3)/2;^ versus x in Fig. [1] also 
demonstrates. A quadratic extrapolation to the physi- 
cal point gives 0.227(38), telling us that SU(3) breaking 
is roughly 20% at the physical point, where x ~ 0.332 
using the PDG values [33 for m+, ra\^ and /^+. We 
compare the result of heavy baryon SU(3) chiral pertur- 
bation theory [1^ for <5su(3) a-s a function of a;, and we 
find that the coefficient of the the linear term in Eq. [2| 
does not vanish. This implies that an accidental cance- 
lation of the low-energy constants is responsible for this 
behavior. 
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0.1 0.2 0.3 0.4 

FIG. 1: (Top) The SU(3) symmetry breaking measure 5su{3)- 
The circles are the measured values at each pion mass, the 
square is the extrapolated value at the physical point, and 
the shaded region is the quadratic extrapolation and its error 
band. 

(Bottom) (5su{3)/2^^ plot. Symbols as above, but the band is 
a constant fit. 



B. Chiral extrapolation 



the values of some low energy constants, the individual 
axial couplings can be fitted with reasonable x^/dof, the 
combined fit seems to fail to describe the data. The con- 
clusion we derive from this is that SU(3) heavy baryon 
chiral perturbation theory fails to describe the lattice 
data. The pion and kaon masses used in our calculation 
arc probably outside the range of validity of this order 
of the pcrturbativc expansion. It is possible that such 
a unified description of the axial couplings requires the 
next order in chiral perturbation theory or smaller pion 
masses. However, since the physical strange quark mass 
is rather large, it seems that SU(3) heavy baryon chiral 
perturbation theory might not be useful in helping ex- 
trapolate the lattice data to the physical point. Here we 
should not forget that one of the reasons for our failure 
to fit with continuum heavy baryon SU(3) chiral pertur- 
bation theory formulas could be discretization errors due 
to taste symmetry breaking in the sea sector. However, 
such errors should be small; they are further suppressed 
to NLO in the expansion since the valence sector domain- 
wall fermions retain chiral symmetry. Hence we do not 
expect that finite lattice spacing corrections will change 
our conclusion. 



Ideally, we should adopt a chiral extrapolation that 
correctly describes the discretization errors due to finite 
lattice spacing as we extrapolate the axial couplings from 
the pion masses we used in our calculation to the physi- 
cal point. However, such a calculation does not exist in 
the literature. Here we consider the next-best available 
option, which is continuum chiral perturbation theory. 

W. Detmold and C.-J. D. Lin worked out the forward 
twist-two matrix element extrapolation at one-loop or- 
der with finite- volume corrections (39j . From their work, 
we see that finite-volume corrections contribute at the 
order of 10^"* or less for the parameters of our calcu- 
lation. Therefore, we will only implement the infinite- 
volume chiral perturbation theory and further simplify 
the formulation by taking mvai = ^sca.- In these formu- 
lae, there are total eight parameters to be determined: 
three SU(3) coupling constants (C, D, F), and five other 
low energy constants (Ac„, Aa„, A/3„, Ajn, Acr„). We 
replace three mass splittings between the decuplet and 
octet with their lattice-measured values and replace / 
and the chiral perturbation theory scale (/x) with the 
pseudoscalar decay constants calculated from the same 
lattice and actions; a similar strategy was adopted by 
NPLQCD and LHPC [2l,[2i,[23,[2i,[2i,[33,[3l|,[32|,[43l. It 
has also been shown that such an approach simplifies the 
mixed action chiral perturbation theory formulas (l8l.[20j. 
We perform a simultaneous fit among multiple axial cou- 
pling constants without making further assumptions. We 
attempted to do simultaneous fitting among all three qa, 
gss and gss, and found the x^/dof for such a fit is of 
the order of 10^. If we limit ourselves to the lightest 
three pion masses, we can reduce x^/dof to the order 
of 10. Although under certain plausible assumptions for 



C. Simple chiral forms 

Since the continuum chiral extrapolation fails to de- 
scribe our data, we take a step back and expand the ax- 
ial couplings in terms of the SU(3) breaking parameter 
X = (m|- — mj)/(47r^/^) as follows: 

gA = D + F + J2C^N^x'' 

71 

n 

5HH = F - + ^ . (3) 

n 

This form reduces to the known SU(3)-symmetric limit 
where the axial couplings are simple linear combina- 
tions of F and D. In addition, not all the C^'' (with 
B G {N, E, S}) are independent parameters. In Sec. lIII Al 
we find that the constraint - 2C^^^ + c4^^ = is 
preferred by the data. This suggests a 4-parameter fit to 
n = 1 order or a 7-parameter fit to n = 2 order if we ig- 
nore possible pion mass dependence of F and D. In order 
to have a reasonable fit form with the smallest number 
of parameters we will focus on the n — 1 case. Figure [5] 
shows our lattice data as a function of (to^/Ztt)^ with the 
corresponding chiral extrapolation; the band shows the 
jackknife uncertainty. 

The x^/dof is 0.83 and the linear fit parameters are 

very poorly determined as C^^ = 0.02(13) and C^^ = 
—0.01(6). This is not surprising, since a small slope is 
seen in all three axial charges. At the physical pion point, 
we find the nucleon axial charge is gA,N = 1-18(4); this 
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is is consistent with what LHPC obtained, 1.23(8), from 
SU(2) chiral perturbation extrapolation [401 • The ex- 
trapolated coupling constants gss and gsH are 0.450(21) 
and —0.277(15) respectively. These numbers are consis- 
tent with the existing predictions from chiral perturba- 
tion theory |4] and large- iVc calculations but have much 
smaller errors. The low-energy chiral parameters arc 
D = 0.715(6) and F — 0.453(5), which arc not consis- 
tent with the recent determination of £> = 0.804(8) and 
F = 0.463(8) using semileptonic decay data and assum- 
ing SU(3) symmetry [ilj. However, the two calculations 
do agree within the range of the SU(3) breaking effect we 
observed in the previous section. 

Since we know that the SU(3) breaking in Eq. [2] is 
quadratic in x, we expect that one needs to go at least to 
71 = 2 in order to capture this effect. Hence, when taking 
the n ~ 2 expansion, it is not surprising to sec the x^/dof 

(ti) 

improves to 0.57. However, the coefficients Cg remain 
poorly determined (most of them are consistent with zero 



within the errorbar). Note that C^^ - 2C^' + C, 
1.9(6) is consistent with what we found for the curvature 
of (5su(3) in Sec. 11101 The final results, D = 0.711(7) 
and F = 0.452(5), are consistent with the n — \ case but 
are better determined. The axial couplings in this case 
are gA = 1.28(6), 5ss = 0.39(6) and gsH = -0.24(4). 
The discrepancy between the results of n = 1 and n = 2, 
which is at one standard deviation, might be taken as an 
indication of the systematic error of such extrapolations. 
However, it is hard to make an honest determination of 
such systematic error without further study at lighter 
pion masses and higher statistics. 

Finally, one can consider that D and F have de- 
pendence as 

D = D„ + CDml/iATi^fl), F = Fo+ cpml/i^n^fl), 

where Dq and Fq are the chiral limit axial couplings. 
Such fit forms have more parameters than we can possi- 
bly determine with our data. To better understand how 

(n) 

strong this dependence is, let us assume m 
Eq. [3] This gives us cd = -0.03(7) and cp = 0.01(5), 
consistent with zero. In both scenarios, the axial cou- 
pling constants are consistent with the n = 1 extrapola- 
tion. The discrepancies in the extrapolated results for all 
fitting forms we used is always at the level of one stan- 
dard deviation. Therefore, we will assign an upper limit 
for the systematic error the same amount as the statis- 
tical error due to the extrapolation and keep the n = 1 
results as our central values. 

The systematic errors due to finite-volume effects are 
expected to be small. LHPC calculated the nuclcon axial 
coupling constant using a chiral extrapolation with finite- 
volume correction j^; less than a 1% effect is observed. 
Furthermore, finite-volume effects (including the E and 
S baryons) are also estimated in Ref. [1^, using heavy 
baryon chiral perturbation theory. The finite-volume ef- 
fects come in at a magnitude no larger than 10~^. Thus, 
we take 1% to be an upper bound for the finite- volume 
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FIG. 2: Lattice data (circles) for qa, gss and gss and chiral 
extrapolation (lines and bands). The square is the extrapo- 
lated value at the physical point. 



systematic error in our calculation. Such effect is negli- 
gible given our statistical errors and systematics due to 
chiral extrapolation. 

The final source of systematic errors is the contin- 
uum extrapolation. Precise estimate of such errors re- 
quires computations at several lattice spacings. How- 
ever, we can estimate the discretization errors which are 
0(a2A|cj3). Taking Aqcd « 300 MeV and the value of 
our lattice spacing, ^ 1588 MeV [l^l, one expects 
to have such a systematic error of the order of 4%. 



IV. CONCLUSION 

In this work, we calculate the axial coupling constants 
for E and S strange baryons using lattice QCD for the 
first time. We do the calculation using 2-|-l-flavor dynam- 
ical configurations with pion mass as light as 350 MeV. 
We discussed various potential chiral extrapolations and 
various sources of systematic errors. We conclude that 
gA = 1.18(4),tat(6)syst, <7EE = 0.450(21),tat(27),y,t and 
gss = -0.277(15)stat(19)syst. In addition, the 5/7(3) 
axial coupling constants are estimated to he D = 
0.715(6)stat(29)syst and F = 0.453(5),tat(19)syst. The ax- 
ial charge coupling of S and S baryons are predicted with 
significantly smaller errors than estimated in the past. 
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